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According to quantum mechanics chiral molecules, that is molecules that rotate the polarization 
of light, should not exist. The simplest molecules which can be chiral have four or more atoms 
with two arrangements of minimal potential energy that are equivalent up to a parity operation. 

Chiral molecules correspond to states localized in one potential energy minimum and can not be 
stationary states of the Schrodinger equation. A possible solution of the paradox can be founded 
on the idea of spontaneous symmetry breaking. This idea was behind work we did previously 
involving a localization phase transition: at low pressure the molecules are delocalized between the 
two minima of the potential energy while at higher pressure they become localized in one minimum 
due to the intermolecular dipole-dipole interactions. Evidence for such a transition is provided by 
measurements of the inversion spectrum of ammonia and deuterated ammonia at different pressures. 

A previously proposed model gives a satisfactory account of the empirical results without free 
parameters. In this paper, we extend this model to gas mixtures. We find that also in these systems 
a phase transition takes place at a critical pressure which depends on the composition of the mixture. 
Moreover, we derive formulas giving the dependence of the inversion frequencies on the pressure. 

These predictions are susceptible to experimental test. 

PACS numbers: 34.10.+X; 03.65.Xp; 34.20.Gj 


I. INTRODUCTION 

The existence of chiral molecules, that is molecules 
that rotate the polarization of light, has so far no univer¬ 
sally accepted explanation. Hund was the first to point 
out in 1927 that such molecules according to quantum 
mechanics should not exist as the corresponding Hamil¬ 
tonian is invariant under mirror reflection [1]. This prob¬ 
lem became known as Hund’s paradox and was discussed 
in the physico-chemical literature for decades. A possi¬ 
ble solution of the paradox can be founded on the idea 
of spontaneous symmetry breaking (SSB) once we take 
into account that the experiments deal with a large num¬ 
ber of molecules and that their interactions are sufficient 
to keep them in the chiral state. This idea was behind 
work done previously [2-8] where the explanation of the 
existence was based on a localization phenomenon. 

Pyramidal molecules, i.e., molecules of the kind XY 3 
such as ammonia NH3, phosphine PH3, arsine ASH3, are 
potentially chiral. Suppose that we replace two of the 
hydrogens with different atoms such as deuterium and 
tritium: we obtain a molecule of the form XYWZ. This 
is called an enantiomer, that is a molecule whose mir¬ 
ror image can not be superimposed to the original one. 
These molecules are optically active in states where the 
atom X is localized on one side of the YWZ plane. In 
the following we shall refer to potentially chiral molecules 
such as ammonia as pre-chiral molecules. Another exam¬ 
ple of pre-chiral molecule, not pyramidal, is deuterated 
disulfane D 2 S 2 (see Fig. 1). 

According to quantum mechanics chiral molecules 
should not exist as stable stationary states. Consider a 
pyramidal molecule XY 3 . The two possible positions of 
the X atom with respect to the plane of the Y atoms are 


separated by a potential barrier and can be connected 
via tunneling. This gives rise to stationary wave func¬ 
tions delocalized over the two minima of the potential 
and of definite parity. In particular, the ground state is 
expected to be even under parity. Tunneling induces a 
doublet structure of the energy levels. 

On the other hand, the existence of chiral molecules 
can be interpreted as a phase transition. In fact, isolated 
molecules do not exist in nature and the effect of the 
environing molecules must be taken into account. This 
interpretation underlies a simple mean-field model that 
we have developed [7] to describe the transition of a gas 
of NH3 molecules from a nonpolar phase to a polar one 
through a localization phenomenon which gives rise to 
the appearance of an electric dipole moment. Even if 
ammonia molecules are only pre-chiral, the mechanism, 
as emphasized in [7], provides the key to understand the 
origin of chirality. 

A quantitative discussion of the collective effects in¬ 
duced by coupling a molecule to the environment con¬ 
stituted by the other molecules of the gas was made in 
[3]. In this work it was shown that, due to the insta¬ 
bility of tunneling under weak perturbations, the order 
of magnitude of the molecular dipole-dipole interaction 
may account for localized ground states. This suggested 
that a transition to localized states should happen when 
the interaction among the molecules is increased. 

Evidence for such a transition was provided by mea¬ 
surements of the dependence of the doublet frequency 
under increasing pressure: the frequency vanishes for a 
critical pressure P cr different for NH3 and ND3. The 
measurements were taken at the end of the 1940s and 
beginning of the 1950s [9-11] but no quantitative theo¬ 
retical explanation was given for 50 years. Our model [7] 
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FIG. 1: (Color online) Molecular structures of ammonia 
(NH 3 , left) and deuterated disulfane (D 2 S 2 , right) in one of 
their two localized states. 


gives a satisfactory account of the empirical results. A 
remarkable feature of the model is that there are no free 
parameters. In particular, it describes quantitatively the 
shift to zero-frequency of the inversion line of NH 3 and 
ND 3 on increasing the pressure. 

In this paper, we extend our model to gas mixtures. 
This case may be of interest, among other things, for 
the interpretation of the astronomical data such as those 
from Galileo spacecraft [12] which measured the absorp¬ 
tion spectrum of NH 3 in the Jovian atmosphere. 

The main results of this work are as follows. 

(i) Also in the case of mixtures there is a quantum 
phase transition between a delocalized (or achiral, or 
nonpolar) phase and a localized (or chiral, or polar) 
phase. At fixed temperature, the crossover between the 
two phases is determined by a properly defined critical 
pressure P CI . 

(ii) We derive formulas expressing the critical pressure 
of the mixture in terms of the doublet splittings of the 
isolated chiral molecules and the fractions of the con¬ 
stituents. For a binary mixture with fractions x\ and xi 
of the two species, we have 


1 
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The paper is structured as follows. In the next section 
we briefly recall the mean-field model for a single gas of 
molecules exhibiting inversion doubling such as ammo¬ 
nia. In particular, we discuss the interaction mechanism 
and the strength of the interaction leading to a localiza¬ 
tion quantum phase transition with the disappearance of 
the inversion line. In Sec. Ill, we describe the well known 
reaction field mechanism which provides exactly the same 
interaction strength. The reason for arguing in terms of 
the reaction field is its simplicity especially in view of its 
extension to gas mixtures. In Sec. IV, we consider mix¬ 
tures of a chiral gas with non polar molecules, such as, for 
instance, NH 3 -H 2 -He. In Sec. V, we consider mixtures 
of chiral gases. For convenience of the reader, we have 
added three appendixes. They give details on the model 
[7] and the comparison with the experiments [9-11]. 


II. CHIRAL MOLECULES AS A CASE OF 
SPONTANEOUS SYMMETRY BREAKING 

We model a gas of pre-chiral, e.g., pyramidal, 
molecules as a set of two-level quantum systems, that 
mimic the inversion degree of freedom of an isolated 
molecule, mutually interacting via the dipole-dipole elec¬ 
tric force. 

The Hamiltonian for the single isolated molecule is as¬ 
sumed of the form — a x AE/2 , where A E is the inversion 
energy splitting measured in a rarefied gas of the consid¬ 
ered species and 


a 


X 


0 1 
1 0 


is the Pauli matrix with symmetric and antisymmetric 
delocalized tunneling eigenstates and ip-: 


Supposing that species 1 is chiral, PcP is the crit¬ 
ical pressure of the pure unmixed species, namely, 
Pc P = A£d/( 27 n), with A E\ being the rarefied-gas 
doublet splitting and 711 given in terms of the tem¬ 
perature and other microscopic parameters, i.e., the 
electric dipole moment and the collision diameter of 
molecules 1 [see Eq. (24)]. If species 2 is non polar, 
then Per = AE , i/( 27 i 2 ) with 712 expressing the polar¬ 
ization of molecules of species 2 by the chiral molecules of 
species 1 [see Eq. (10)]. If species 2 is chiral, then Per = 
AE 2 /( 2722 ) with 722 expressing the self-consistent inter¬ 
action among the molecules of species 2 [see Eq. (24)]. In 
both cases, the inverse critical pressure of the mixture is 
the fraction-weighted average of the inverse critical pres¬ 
sures of its components. The result is easily generalized 
to mixtures made up of an arbitrary number of compo¬ 
nents. 

(iii) We derive formulas giving the dependence of the 
inversion frequencies on the pressure [see Eqs. (12) and 
Sec. VB], These results are susceptible to experimental 
test. 


V2 (O’ ^ V2 (- 1 )' 

Since the rotational degrees of freedom of an isolated 
molecule are faster than the inversion ones, on the time 
scales of the inversion dynamics, namely, h/AE , the 
molecules in a gas feel an attraction arising from the angle 
averaging of the effective dipole-dipole interaction at the 
actual temperature of the gas [13]. The localizing effect 
of the dipole-dipole interaction between two molecules i 
and j can be represented by an interaction term of the 
form —Qijcrf (g> crj, with g^ > 0 , where 



is the Pauli matrix with left and right localized eigen¬ 
states ipL and tpR 
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The Hamiltonian for N interacting molecules then reads 
as 


AF N N N 

H = - 2 E<-E E a) 

i=l i=l j=*+l 

For a gas of moderate density, we approximate the 
behavior of the N 1 molecules with the mean-held 
Hamiltonian 


AF G 

h(iP) = -—<J*--(iP,cj*iP)c z , (2) 

where ip is the mean-held molecular state to be deter¬ 
mined self-consistently by solving the nonlinear eigen¬ 
value problem associated with Eq. (2) and the normaliza¬ 
tion condition (ip, ip) = N. The scalar product between 
two Pauli spinors ip and <p is defined in terms of their two 
components in the standard way 

(lp,(p) = 1p l( pi + 1p 2 <p2- 

The parameter G accounts for the effective dipole in¬ 
teraction energy of a single molecule with the rest of the 
gas. It must be identihed with a sum over all possible 
molecular distances and all possible dipole orientations 
calculated assuming that, concerning the translational, 
vibrational, and rotational degrees of freedom, the N 
molecules behave as an ideal gas at thermal equilibrium 
at temperature T. This assumption relies on a sharp sep¬ 
aration (decoupling) between these degrees of freedom 
and the inversion motion [14]. 

We now discuss the calculation of G following Keesom 
[13]. An alternative derivation of G will be shown later 
based on the so called reaction held mechanism. 

Let us dehne an effective interaction potential V(r) 
between two electric dipoles at distance r via the formula 

e -V(r)/k B T = / e -V(r,n u u. 2 )/k B Ty 


where 


V(r,n i,/i 2 ) 


Mi • M 2 ~ 3(/Xj ■ r)(r ■ fi 2 )r 2 
47T£of 3 


is the standard dipole-dipole interaction potential and 
(...) denotes the average over all possible orientations 
of the vectors /ii and M 2 - Assuming that V(r) <C fegT 
and observing that (V(r,[M 1 , 1 * 2 )) = 0, we approximately 
evaluate the effective potential as 


2 2 

= mim 2 
3(47T£ 0 ) 2 fcsTr 6 ’ 

In the case of a gas of molecules with dipole moments p, 
the effective interaction energy of one molecule with the 


TABLE I: Rarefied-gas energy splitting A E, electric dipole 
moment p, and collision diameter d measured for different 
pre-chiral molecules [14, 15]. In the fourth column, we report 
the critical pressure P cr evaluated by Eq. (4) at T = 300 K. 



A E (cm -1 ) 

M (Debye) 

d( A) 

Per (atm) 

nh 3 

0.81 

1.47 

4.32 

1.69 

nd 3 

0.053 

1.47 

4.32 

0.11 

d 2 s 2 

10 -9 

1.56 

5.97 

4.3 x 10 -9 


rest of the gas amounts to —G, where 

r°° „ 4 

G= / p(r)—--rw;——-^-47rr 2 dr 

J d ^ ; 3(4^£ 0 ) 2 fc B Tr 6 

_ 4tt p 4 P 

9 (47r£ofcsT) 2 (i 3 


(3) 


To get the above result we used p(r) = P/ksT , the ideal 
gas density at pressure P and temperature T. Note that 
d is the minimal distance at which two molecules can be 
found, namely, the so called molecular collision diameter. 
At fixed temperature, the effective interaction constant 
G increases linearly with the gas pressure P. 

The nonlinear eigenvalue problem associated with 
Eq. (2) is solved in detail in Appendix A. There exists 
a critical value of the interaction strength, G cr = AE/2, 
such that for G < G cr the mean-field eigenstate with min¬ 
imal energy is ip = y/Nip+, with delocalized molecules, 
whereas for G > G cr there are two degenerate states of 
minimal energy. These are chiral states transformed into 
each other by the parity operator a x approaching the lo¬ 
calized states ipL, ifR for G G cr . We thus have a quan¬ 
tum phase transition between a delocalized (or achiral, 
or nonpolar) phase and a localized (or chiral, or polar) 
phase. In view of the dependence of G on P, we can 
define a critical pressure at which the phase transition 
takes place 


9 AEd 3 (47T£ 0 k B T) 2 

cr 87r p 4 ‘ U 

In Table I we report the values of P cr calculated for dif¬ 
ferent pre-chiral molecules. 

In Appendix B we show that in the delocalized 
phase the inversion angular frequency of the interacting 
molecules depends on the pressure as 


^ = AE \j 1 - jr- ( 5 ) 

This formula is interesting as it expresses the ratio of 
two microscopic quantities, Fuj and A E, as a universal 
function of the ratio of the macroscopic variables P and 
P cr . Furthermore, it provides a very good representation 
of some spectroscopic data, namely, the shift to zero fre¬ 
quency of the inversion line of NH 3 or ND 3 , as it appears 
from Fig. 11, see Appendix C for more details. 
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III. THE REACTION FIELD MECHANISM 

The mean-field coupling constant G introduced in Sec¬ 
tion II can be estimated, obtaining the same result, us¬ 
ing a different argument which can be easily general¬ 
ized. This is the reaction field mechanism widely used 
in physics and chemistry [16, 17]. 

Let us consider a spherical cavity of radius a in a ho¬ 
mogeneous dielectric medium characterized by a relative 
dielectric constant e r . An electric dipole /z placed at the 
center of the cavity polarizes the dielectric medium in¬ 
ducing inside the sphere a reaction field R proportional 
to fj,: 


r _ 2(g r 1) M 

2e r + 1 47 T£ 0 a 3 ' 

As a result, the dipole acquires an energy 




R 


£ r — 1 / 1 2 

2e r + 1 47T£ 0 a 3 ’ 


By using the Clausius-Mossotti relation 


e r — 1 
£r T 2 


L , 

p{a + i 


( 6 ) 


where a is the molecular polarizability and a D = 
^ 2 /(3 £ o ksT) is the Debye (orientation) polarizability, 
and observing that for a chiral gas a D a (for in¬ 
stance, in the case of NH 3 we have a ~ 2 A whereas 
a D ~ 217 A 3 at T = 300 K), we get 

„ _ 47r p 4 P 

9 ( 4:7T£ 0 kBT) 2 a 3 


In all the above expressions we have £ r ~ 1 so that we 
approximated £ r + 2 ~ 2 £ r + 1 ~ 3. 

We exactly have § = — G provided a = d, where 
d is the molecular collision diameter introduced in 
Eq. (3) as the minimum distance between two interact¬ 
ing molecules. In fact, the size of the spherical cavity is 
not arbitrary. Microscopic arguments [18, 19] show that 
a must be identified with the radius of an effective hard 
core, namely, the molecular collision diameter d. 


IV. MIXTURES OF A CHIRAL GAS WITH 
NON POLAR MOLECULES 

By using the reaction field approach to the coupling 
constant G, we now generalize the mean-field model of 
Section II to the case of a mixture of chiral with non polar 
molecules. Consider a gas mixture of two species labeled 
1 and 2. In this case, the Clausius-Mossotti relation reads 
as 


£ r — 1 
£ r ~t~ 2 


1 

3 


(pi (oq + a?) + p -2 (a 2 + a!?)) ■ 


If species 1 is pre-chiral, i.e., polar, oq can be neglected 
with respect to af and if species 2 is non polar a§ = 0 . 
The energy acquired by one molecule of the polar species, 
having electric dipole moment /zi, then is 


1 


Mi 


£\ — — 7 : Pi oP —-^ 5 — 

3' 1 4:'K£ 0 dl 1 


1 


Mi 

3 /32Q2 47 r g 0 di2 ’ 


(7) 


where dn and g?i 2 are the diameters of the molecular 
collisions 1-1 and 1 - 2 , namely, dn = 2 ri and dr 2 = rq + 
r' 2 , rq and 7*2 being the hard sphere radii of the molecules 
1 and 2. If the density pi is such that <C P 2 & 2 , the 
above formula reduces to 

* _ 1 Mi 

^~~3 P2a2 4^dl 2 - 


This formula approximates the case of a single polar 
molecule immersed in a non polar gas. 

For a mixture of chiral and non polar molecules, the 
mean-field molecular state of the chiral species 1, 7 /q, is 
determined similarly to the case of a single chiral gas. 
The only degree of freedom of the non polar molecules 
is the deformation which, in turn, is proportional to the 
electric dipole moment of the chiral molecules. Therefore, 
we assume the following nonlinear Hamiltonian similar to 
Eq. (2) 

hi(ipi) = -(8) 

The mean-field molecular state 7 /q is normalized to the 
number of molecules of the species 1, namely, (V’i, V’i) = 
Ni. The parameter Gi accounts for the effective dipole 
interaction energy of a single chiral molecule with the 
rest of the gas constituted by polar and non polar species. 
According to the reaction field arguments, we have just 
to identify G\ = —£\, where £\ is given by Eq. (7). Thus 
we see that Gi is the sum of two contributions 


G\ = Gn + G 12 , 


where 


G u = 


47T 


pfpl 


9 (ATT£okBT) 2 d\ 


= I 11 P 1 


(9) 


11 


is the energy provided by the interaction with the chiral 
species and 


= 1 02 PIP2 

12 3 4:7T£okBTdi2 


712-F2 


( 10 ) 


the energy due to the presence of the non polar molecules. 
In the above expressions we used the ideal gas relations 
Pi = Pi/ksT and p 2 = JAAbT" , Pi and P 2 being the 
partial pressures of the species 1 and 2 . 

The analysis of the mean-field Hamiltonian ( 8 ) is iden¬ 
tical to the case of a single chiral gas. We have a lo¬ 
calization phase transition when G 1 becomes equal to 
Gi cr = AEi/2. The transition can be considered as a 
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FIG. 2: (Color online) Critical pressure P cr as a function of 
the chiral species fraction x for NH^-He (bottom) and D 2 S 2 - 
He (top) mixtures at T = 300 K. Note that pressure units are 
atm in the case of NH 3 He and natm in the case of D 2 S 2 -He. 


function of the total pressure P = Pi+ P 2 oi the mixture 
and of the fractions of the two species x\ = P\/P and 
x 2 = P 2 / P. In this case, an important prediction is that, 
instead of a unique critical pressure, we have a critical 
line parametrized by X\ (or x 2 = 1 — Xi): 

Pcr{Xl) = 2cr l7ll +2(l 1 -* 1 )7i 2 ’ (11) 

where 7 n = Gn/Pi and 7 i 2 = G\ 2 /P 2 with Gn and 
G 12 as in Eqs. (9) and (10). 

The inversion frequency of the chiral species 1 can be 
evaluated as in the case of a single chiral gas. The result 
is formally identical to that in Eq. (5), 


Kuj 1 — AEi 



P 

-PcrOEl)’ 


( 12 ) 


with the angular frequency uq which is now a function 
of the pressure P of the mixture and of the fraction x\ 
of the chiral species. However, for any fixed value of aq, 
the ratio fiuJi/AEi is still a universal function of P/P cr 
independent of the nature of species 1 . 

In Fig. 2, we show the variation of the critical pressure 
as a function of the fraction x of the chiral species for 
the mixtures NH 3 -He and D 2 S 2 - He. The value of P cr 
at x = 1 coincides with that reported in Table I for a 
pure chiral gas, whereas for x —> 0 the critical pressure 
approaches a much higher value. The ratio 

lim x -» 0 PcrO) _ 4 tt ^l/^endh) d\ 2 
Pcr(l) 3 k B T a 2 

amounts to about 400 in the examples of Fig. (2). 

The above analysis is immediately generalized to a 
mixture containing one chiral gas and an arbitrary num¬ 
ber of non polar gases. As an example, for a ternary 


mixture (as before, species 1 is the chiral one) we have a 
surface of critical pressures 

P CT (xi,x 2 ) 2 a; l7ll + 2a; 27 i2 + 2(1 - x\ - 2 : 2 ) 713 ’ 

where 713 = G 13 /P 3 with G 13 given by an expression 
analogous to Eq. (10). This term expresses the polariza¬ 
tion of the molecules of species 3 by the chiral molecules 
of species 1 . 

The above results may be of interest, among other 
things, in the interpretation of astronomical data such 
as those from Galileo spacecraft [12] which measured the 
absorption spectrum of NH 3 in the Jovian atmosphere. 
Note that, since the atmosphere of Jupiter is mostly made 
of molecular hydrogen and helium in roughly solar pro¬ 
portions whereas ammonia has a tiny fraction, the criti¬ 
cal pressure of the mixture is drastically different if com¬ 
pared with the 1.69 atm value of ammonia. For fractions 
2 ;nh 3 — 0.01, Xh 2 — 0.86 and x Be ~ 0.13, we estimate 
P cr ranging from 15 atm at T = 100 K to 165 atm at 
T = 400 K. This implies that no ammonia inversion line 
shift can be observed in microwave absorption measure¬ 
ments performed in the atmosphere of Jupiter. 


V. MIXTURES OF CHIRAL GASES 

It is possible to extend our mean-field approach to gas 
mixtures made of several chiral species. We start consid¬ 
ering the simplest case of two species. 

In Sec. V A, we illustrate the model and evaluate the 
stationary molecular states of the mixture. In a first 
reading, one can skip the discussion after Eq. (18) and 
look at the results summarized in Table II. As in the 
case of a single chiral gas, the mixture undergoes a lo¬ 
calization phase transition at a critical pressure P cr . For 
0 < P < P C r, the lowest-energy molecular state of the 
mixture, named d+y, corresponds to molecules of both 
species in a delocalized symmetric configuration. For 
P > P cr , new minimal-energy molecular states appear 
with twofold degeneracy. These states, named and 
crr, correspond to molecules of both species in a chiral 
configuration of type L or R. Equation (25) provides a 
simple formula for P cr . The inverse critical pressure of 
the mixture is the fraction-weighted average of the in¬ 
verse critical pressures of its components. 

In Sec. VB, we evaluate the excitation energies of a 
two-species mixture by using the linear-response theory. 
The results are summarized by Eq. (36). We have two 
normal modes which at low pressure reduce to the inver¬ 
sion of the component species. By increasing the pres¬ 
sure, the energy of both modes decreases, the smaller 
one vanishing at P = P cr , the larger one approaching for 
P P cr the asymptotic value given by Eq. (37). An 
example of this behavior is shown in Figs. 5 and 6 in the 
case of a NH 3 -ND 3 mixture. 

Finally, in Sec. V C, we discuss the consistency of the 
developed theory by considering the absorption coeffi- 
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cient and two limits in which it must reduce to that of a 
single chiral gas, namely, the limit of zero fraction of one 
species and the limit of equal species. We also discuss 
the possibility of experimental tests for the considered 
NH 3 -ND 3 mixture. 

All the above results can be generalized to mixtures of 
several chiral gases. 

A. Molecular states 

For a mixture of two chiral gases, the Clausius- 
Mossotti relation reads as 

£r 1 1 / 73 

c _l_ 2 = 3 w lQ b + P 2a 2 ) > 

where the Debye polarizations af = /zf/(3£o k B T), 1 = 
1 , 2 , are given in terms of the molecular electric-dipole 
moments, pi and p , 2 , of the two species. According 
to Eq. ( 6 ), due to mutual interactions the molecules of 
species 1 and 2 acquire the energies 

_ 47t /xf/xfPaq 4 tt plp 2 2 Px 2 

1 9 (47T£o kBT) 2 di 1 9 (47T£q ksT^d^ ’ 


2, separately. In fact, the molecules of the two species 
can not be transformed into each another. 

As discussed in the case of a single chiral gas (see Ap¬ 
pendix A) the molecules of the species 1 and 2 will occupy 
the eigenstate having the lowest energy 

2 A /■ 

S{ipi,ip2) = ~Y^ 

i=1 
2 2 

“EE 2N~.^ u o'lV’iXV’j, Vjipj)- (16) 

i =1 j =1 1 

In virtue of the symmetry property G 12 /N 2 = G 2 \/Ni, 
the above energy is a constant of motion of the infinite¬ 
dimensional Hamiltonian system 

l^Vd = /q(tq,t’2)t-i. 

Notice, in fact, that <5£('i/’i, )ipi- 

To solve the eigenvalue problem (15), we write 


47T /i 2 /-t|Px 2 47T /if/zfPa; 1 

9 {A'K£Qk B T) 2 d\ 2 9 (4jr£ok B T) 2 d2 2 


As usual, we used the ideal gas relations pi = Pi/k B T, 
1 = 1 , 2 , and introduced the fractions X\ = P\/P and 
X 2 = P 2 /P 1 where P is the total pressure of the mixture. 

We describe the inversion degrees of freedom of the two 
species by mean-field molecular states , ip 2 normalized 
to the number of molecules of the corresponding species, 
(it>i,4>i) = Ni, (%/j 2 ,^ 2 ) = N 2 . Here, Vh and ^2 are Pauli 
spinors and the scalar product is defined as in Section II. 
By writing <£) = — JA =1 G lt j, i = 1,2, where 


4?t tftfjPxj 

9 (4:TT£ 0 k B T) 2 d^ 


(13) 


and reasoning as in Sec. II, we introduce the following 
set of two nonlinear Hamiltonians: 


hi(ipi,ip 2 ) 


h 2 {tp\, ^ 2 ) 


APi 

-< 

2 

AE 2 


E 

t=i 


-E 


^7 


(14a) 

(14b) 


Pauli operators now have a label 1 = 1,2 relative to the 
species they refer to. Stationary states are obtained by 
solving the eigenvalue problem 


hi{ipi,ip 2 )ipi = AiV’i, (15a) 

h 2 {’>pi,^2)i’2 = A 2 V’ 2 , (15b) 

where A 1; A 2 are the Lagrange multipliers associated with 
the conservation of the particle number of species 1 and 


ipi = 1 = 1 , 2 . (18) 

The coefficients at and 5,; can be chosen real and are 
normalized to the number of molecules of each species, 
af + b 2 = Ni. Up to an irrelevant global sign, we can 
also choose ai and 02 non-negative. By inserting the 
expressions (18) into Eq. (15), we find that ai,bi,a 2 ,b 2 
are the solutions of the system of equations 

2 

A Eicubi = ^2 2 g i:j (af - 6 f) ajbj, 1 = 1,2, (19) 

i =1 


where g t j = Gij/Nj. Once a solution of this system is 
found, the corresponding energy can be evaluated as 

2 AF. 22 

£ = - E (°< ~ 6 <) ~ E E 2gij a ibidjbj. 

i—1 i—1 j—1 

It is convenient to discuss the solutions of Eq. (19) in 
function of the parameter P £ [0, 00 ). In fact, the matrix 
elements g t j are proportional to the pressure. For any 
value of P , Eq. (19) admits solutions such that ai&i = 0 
and a 2&2 = 0. There are four different solutions satisfying 
this condition and the normalization rule (see Table II). 
The solution 

a\ = \/ bi = 0 , a 2 = \J N 2 , b 2 = 0 ( 20 ) 

is that with the lowest energy £ = —N(x\AEi/2 + 
X 2 AE 2 / 2 ). We name this solution d ++ because it cor¬ 
responds to molecules of both species in a delocalized 
symmetric configuration. 
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TABLE II: Stationary molecular states for a mixture of two chiral gases with fractions xi = N\ /N and X 2 = N 2 /N. N = N 1 +N 2 
being the total number of molecules. The states are classified in column 1 as delocalized ( d ) even (+) or odd (—) and chiral 
(c) left (L) or right ( R ). This classification results from the decomposition of the two-species molecular states in terms of the 
delocalized eigenfunctions of a x via the coefficients ai, 61 , a 2 , 62 as in Eq. (18). The values of ai, 61 , 02 , 62 are found by solving 
Eq. (19) with the normalization conditions af + bf = Ni and a 2 + b 2 = N 2 - The chiral solutions exist only for P > P CI and 
are defined in terms of two parameters qi and 52 shown in Fig. 4. The energy per unit molecule associated with each state [see 
Eq. (16)] is reported in the last column. For P < P CI the solution of minimal energy is d++ while for P > P CI there are two 
degenerate solutions of minimal energy, namely, cll and crr. Other higher-energy states bifurcating at pressures larger than 
Per, see Fig. 3, are not reported in this table. 



ai 

bi 

0,2 

62 

£/N 

d++ 

Vm 

0 

VW 2 

0 

x\AE\ X2AE2 

2 2 

d+- 


0 

0 

^/W 2 

xiAEi | X2AE2 

2 1 2 

d- + 

0 

VNi 


0 

1 x\AE\ X2AE2 

'2 2 

d__ 

0 

■v/W 

0 

■s/W 2 

1 xiAEi | X2AE2 
' 2 ' 2 

cll \ 

Jit 0 - + qi) \J 


\J ^(1 + 52) 

Vf a-®) 

a 771 Xj Xi \ /l — < 7 ? . /l — n P 

xiqiAEi * JV V ' 3 

2 -^i 2 Z— /i,j 2 

crr \ 

i + 91) 

/^( l- 9 i) 

\J ^(1 + 52) 

-\/aM 1-92) 

a j—i Xa Xa \/l — Q? a / l — 4 P 

XiqiAEi 1 J V V ,%3 

2 -^i 2 2 

If CI161 

= 0 and a 2 b 2 7^ 0 

or, vice versa, 

0161 7! 0 and 

where we set 



a 2 b -2 = 0, Eq. (19) has no solutions. In fact, if ai 6 i = 0 
then a\ — b\ = ±N\ and Eq. (19) for i = 1 gives 0 = 
±2gi 2 Nia 2 b 2 which can not be solved for a 2 b 2 7 I 0 . 

Equation (19) admits solutions with ai&i 7 I 0 and 
a 2 b 2 7 I 0 if and only if 

[AEi — 2gn(af — bf)] [AE 2 — 2g 22 (a 2 — b 2 )] 

- ^912921 (a? - 6?)( a 2 - l> 2 ) = 0 . ( 21 ) 

In fact, think of Eq. (19) as a linear homogeneous system 
of equations in the two unknowns 01&1 and a 2 b 2 . This 
system admits a nontrivial solution (aib 1 , a 2 b 2 ) 7 ^ ( 0 , 0 ) 
if and only if its determinant is 0, namely, Eq. (21). Any 
non trivial solution is necessarily of the form (1161 7 ^ 0 
and a 2 b 2 7 ^ 0 because we have already established that 
Eq. (19) does not have solutions with just one term be¬ 
tween aifci and a 2 b 2 equal to 0 . 

If ( 01 , 61 , 02 , 62 ) is a solution of Eq. (19), that is true 
also for ( 01 ,- 61 , 02 ,- 62 ) and these two solutions have 
the same energy. Therefore, each eigenstate obtained for 
0161 7 ^ 0 and a 2 b 2 ^ 0 has a twofold degeneracy. 

Observing that 511522 — 512521 = 0, Eq. (21) can be 
rewritten as 

45ii a\ AE 2 + Ag 22 a 2 AE 1 


N 

2/ij — -p9ij — 


N ^ 47t 

N^ &ij = 


9 2 


9 (47T £Qk B T) 2 d% 


(24) 


For x\ = 1 or x 2 = 1, Eq. (23) reduces to the critical 
pressure (4) for the single species 1 or 2. 

In conclusion, similarly to the case of a single chiral 
gas, it exists a critical pressure P cr and a quantum phase 
transition taking place at P = P cr . For 0 < P < P cr , the 
lowest-energy stationary state of the mixture is the delo¬ 
calized d ++ configuration given by Eq. (20). At P = P cr 
a bifurcation occurs and two new energy-degenerate sta¬ 
tionary states appear (see Fig. 3 for an example). We call 
them cll and Crr because they correspond to molecules 
of both species in a chiral configuration of type L or 
R (see Table II), achieving a complete localization for 
P ^ Per- The value of P cr is given by Eq. (23). This 
formula says that the inverse critical pressure of the mix¬ 
ture is the fraction-weighted average of the inverse critical 
pressures of its components, 


1 

Per 


= £< 

i=1 


1 

p(f 

1 cr 


_p(*) = 


A Ej 

2 ^a 


(25) 


— 2g\iN\AE 2 + 2g 22 N 2 AE\ + AE\AE 2 . (22) 

The left-hand-side of Eq. (22) is manifestly positive, 
therefore we must conclude that also the right-hand-side 
is so. In terms of P, remember that 5,7 oc P, this condi¬ 
tion implies that the pressure has to be larger than the 
critical value 


The difference between the lowest-energy stationary 
solutions of the two phases, namely, d++ and Cll,crr, 
can be entirely encoded in the parameters 


/ 1 , P < P cr , 

\ 2af /Ni - 1, P > P cr , 


* = 1 , 2 , 


Pc, 



C 2‘Kplii 

AEi 


-1 


(see Table II). An example of the behavior of q± and q- 2 
(23) as a function of the pressure is shown in Fig. 4. Each 
parameter qi vanishes for P —> 00 . More precisely, by 



: t s 




1 


1.0 
0.8 h 
0.6 

0.4 

0.2 

0.0 

I 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 


1 


o 


CLL c rr 


3 

P/Pcr 


1_1_1_1_1_1 J 

0 1 

L 1_1_1_1_1_1_1_1_1 

2 3 

P/Pcr 

1 1 1 1 

4 

i 1 _ 1 _ 1 

5 

i 1 _ 1 

6 

: 1 

' . , C LL Crr 




■ VT 






FIG. 3: (Color online) Coefficients a\/\fN\ and 02 /VIV 2 
defining the molecular states tpi , ■02 of a two-species mixture 
as a function of P/P cr . The coefficients have been obtained 
by solving numerically Eq. (19) in the case of a NH3-ND3 
mixture (ND3 is species i = 2) with X 2 = 0.1 at temperature 
T = 300 K. The coefficients bifurcating at P = P CT and tend¬ 
ing for P P cr to 1/\/2 (dashed line) describe two degenerate 
molecular states, named cll and crr, with molecules of both 
species in a chiral configuration of type L or R. Notice fur¬ 
ther bifurcations appearing at P ~ 4P cr , they correspond to 
stationary states of higher energy. 


using Eq. (19) it is simple to show that at the leading 
order for P large one has 
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P/Pcr 

FIG. 4: (Color online) Parameters qi,q 2 encoding the differ¬ 
ence between the delocalized symmetric state d++ and the 
chiral states cll and crr of Table II evaluated numerically 
for the same mixture of Fig. 3. 


ing an electro-magnetic radiation. In detail, we mod¬ 
ify the Hamiltonians (14) by adding for each species 
* = 1,2 a perturbation of the form /qF(i)of with the 
field F(t) = F cos (uit) sufficiently small so that effects of 
order 0{F 2 ) are negligible. The perturbation induces a 
modified time evolution of the unperturbed ground state 
of the system which we evaluate as follows. 

For P < P cr , the ground state of the mixture at zero 
external field is the symmetric delocalized state d++ de¬ 
fined by the coefficients ( 01 , 61 , 02 , 62 ) given by Eq. (20). 
For P > P CT , the unperturbed mixture has two lowest- 
energy degenerate stationary states, namely, the states 
cll and crr, defined by the two sets of coefficients 
( 01 , 61 , 02 , 62 ) and ( 01 ,- 61 , 02 ,- 62 ) given in Table II in 
terms of the parameters < 71 , < 72 - For a racemic mixture, 
which is usually the case of concern, we may assume that 
the unperturbed ground state of species * = 1,2 in the 
localized phase is represented by the stationary density 
matrix obtained by summing with equal weights the pure 
chiral states of types L and R 


P 


(0 


Qi = 


Sj = l 7 ij x j 


p 


(26) 


(i) 

where P cr is the critical pressure of the zth species as 
given by Eq. (25). 


B. Excitation energies: Normal modes, inversion 
frequencies 

As in the case of a single chiral species, we can de¬ 
termine the excitation energies from the ground state 
of a binary mixture by evaluating the linear response 
to an external time-dependent perturbation represent- 


„( 0 ) 


2 ( 6- ) ( CH h ) + 2 ( -bi ) ( Ui bl ) 

|fili + Qiof). (27) 


Since for P < P cr we have q-\ = q -2 = 1, Eq. (27) rep¬ 
resents the ground-state density matrix of the species 
i = 1,2 also in the delocalized phase. In conclusion, we 
can analyze simultaneously both phases P ^ P cr assum¬ 
ing that the unperturbed lowest-energy stationary state 
is given by the density matrices of Eq. (27). 

Under the effect of the perturbation, the time evolution 
of the density matrices of the two species is governed by 
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the equations 


where 


ih^Biit) = [hi(Qi(t), Q2(t)) + HiF(t)af, £>i(t)], (28a) 

ih^g 2 (t) = [h 2 {gi{t), g 2 {t)) + H 2 F{t)a%, g 2 (t)}, (28b) 

where the mean-field Hamiltonians hi(gi,g 2 ), i = 1,2, 
are the same as Eq. (14) only rewritten in terms of the 
density matrices g\ , g 2 


C(u) — £i/x 2 [h 2 0J 2 — AFf + £ 2 (722 — 712 ) P\ 

+ [ftV - A E\ + a ( 7ll - 721 ) P\ , (33) 
D(u) = ( h 2 u 2 - AF 2 + & 711 -P) 

X [h 2 UJ 2 — AE% + ^2722 P) — a^27l2721 P 2 1 

(34) 

a = 2 XiQiAEi, i = 1,2. (35) 


Mfl, £>2) =- 9i i tr ( a j Q i) a i ■ ( 29 ) 

1=1 

The density matrix of the ith species at time t is written 
as 

Qi(t) = q[ 0) + Xi(t)a f + yi(t)ay + Zi(t)a f. (30) 

Neglecting second order terms in the perturbation, the 
Bloch vector (xi{t),yi(t), Zi{t)) is determined by the sys¬ 
tem of differential equations 

hx' z (t) = 0 

hy[{t) = AEiZi(t) - Y?j=i 2N i q i g ij z j {t) + Wq, ; /qF(f) , 
hz'iit) = -A Eiyi{t) 


The poles of the linear-response function (32) provide 
the excitation energies from the ground state predicted 
by our mean-field model [21]. It is worth to explore first 
two limit cases. 

For noninteracting molecules, namely, 7 ^ = 0, i. j = 
1 , 2 , we have 


xM 


P_ ( iuA 
e 0 \h 2 u 2 - AE 2 


H 2 UJ 2 - AE 2 ) ' 


As expected, the response to an electromagnetic radia¬ 
tion incident on the mixture would, in this case, provide 
excitation energies AE\ and AE 2 corresponding to the 
inversion of isolated molecules of species 1 and 2 . 

For a gas with just one chiral species, let us say x 2 = 0, 
Eq. (32) reduces to 


obtained by plugging the definition (30) into Eq. (28). 

We are interested in evaluating the electric dipole mo¬ 
ments of each species i = 1 , 2 , namely, 

Pi(t) = yi tr (a- gi(t)) = 2/qZi(f). 

From the differential equations for the Bloch vector, we 
find that the Fourier transforms pi(u}),p 2 (ijj) are given 
by the system of algebraic equations 

2 

(h 2 u) 2 - A Ef) p^tx) + 2N i q i AE i g ij p j (u) 

3 =1 

= 2p, 2 N i q i AE i F(u), i = 1,2, (31) 


where F(w) is the Fourier transform of the applied field. 
The frequency-dependent electric susceptibility of the 
?’th species is finally evaluated as the ratio between the 
Fourier transform of the polarization density of that 
species and the Fourier transform of the applied field, 


XiM 


Pi(u)/{N/p) 

£ 0 F(uj) 


where JV/p is the volume occupied by the mixture and 
£0 the vacuum dielectric constant [20]. As usual, the 
density of the mixture at pressure P and temperature T 
is assumed to be p = P/ksT. By determining pi(w) and 
p 2 {pj) from Eq. (31), we find that the susceptibility of the 
mixture, % = \i + X2> is given by 


xM 


P C(co) 
£0 D{w) ’ 


(32) 


/ \ = PP\ _ 2qiA.Fi _ 

X[UJ> £0 h 2 L0 2 - (AE 2 -2q 1 AE 1111 P)' 

In this case we also have [22] qi = 1 for P < P cr and q\ = 
AFi/(2Gn) for P > P cr , where Gu = 711 P, so that we 
recover the result of Eq. (Bl). The same conclusion is 
obtained if we consider a single chiral gas as the mixture 
of two identical species with arbitrary fractions X\ and 
x 2 such that x\ + x 2 = 1 . 

For a general mixture, no partial cancellation of zeros 
between numerator and denominator of x(w) takes place 
as in the above two limit cases. The poles of the electric 
susceptibility coincide with the zeros of P(w). By using 
the property 711722 — 712721 = 0 , we find that the bi¬ 
quadratic equation D(oj) = 0 has solutions with modulus 
uJ± given by 


2H 2 uJ 2 ± =A 12 ± \Ja 12 2 -AB 12 , (36) 

where 

A i2 = AF 2 + AE% — (£i7ii + £2722 )P, 

B\ 2 = AF 2 AFf — (£1711 AFf + £2722 AF 2 )P. 

Clearly, the excitation energies of the mixture hiJ± are 
functions of the pressure. 

At P = 0, the normal modes represent the inversion of 
isolated molecules of species 1 or 2 , namely, 

tuA±{ 0) = A F i± , 

where i+ = 1 and i_ = 2 if AFi > AF 2 and i+ = 2 and 
= 1 if AE 2 > AFi. 
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In the delocalized phase 0 < P < P cr , we have 

B 12 (P) = AE*AE% (1 - P/P cr ), 

in virtue of the fact that q\ = q 2 = 1. For the same 
reason we also have that and £2 are constant in this 
phase so that A\ 2 (P) decreases linearly with P. As a 
result, we find that hco±(P) decrease from their P = 0 
value according to 


h*j±(P) ~ AE i± ^l ^ (i±) • 

These approximated relations, which are merely the in¬ 
version frequency laws for the pure species [see Eq. (Cl)], 
are valid for P small. The larger is the fraction of the 
species i± , the wider is the range of P where the approx¬ 
imation is adequate. 

At P = P cr , we have Bi 2 (P cr ) = 0 and, according to 
Eq. (36), the energy vanishes with an infinite slope. 
The energy of the other mode attains the finite value 
hjJ + = yj A 12 (Per) ■ 

In the localized phase P > P cr , where % = 2af/Ni — 1, 
i = 1,2, with d\ and a\ satisfying Eq. (22), we have 
B\ 2 {P) = 0 for any value of P. It follows that = 0 
throughout this phase. The energy of the other mode 
persists in decreasing and for P P cr approaches the 
asymptotic value 


fioJ™ = 


APi7l2^2 


A 2 ?f 721*1 


7ii*i + 712*2 722*2 + 721*1 


(37) 


This value corresponds to 7^12 evaluated with the lead¬ 
ing order expressions of qq, q 2 as in Eq. (26). 

In Figs. 5 and 6 we show the excitation energies fuD± 
calculated for a NH 3 -ND 3 mixture with 10% and 0.01% 
deuterated ammonia fraction, respectively. 


C. Consistency of the model: Absorption 
coefficient 

We have already discussed the consistency of our model 
by considering two limiting cases in which it must recover 
the single gas case. In fact, we have observed that the 
susceptibility of a binary mixture reduces to that of a sin¬ 
gle chiral gas in the limit of fraction of one species tending 
to zero or in the limit in which the parameters defining 
the two species become equal. Here, we study in deeper 
details this limiting behavior by taking into account the 
absorption coefficient associated with the evaluated sus¬ 
ceptibility. 

According to the Beer-Lambert law, the intensity of a 
monochromatic radiation of angular frequency u> propa¬ 
gating for a length t into the mixture is decreased by a 
factor d ue to the excitation of the normal modes. 



P (atm) 


FIG. 5: (Color online) Energies tuJ5± of the excitation normal 
modes as a function of pressure P in the case of a NH 3 - 
ND 3 mixture (ND 3 is species i = 2) with x 2 = 0.1 at T = 
300 K. The horizontal dotted line is the asymptotic value 

huj+. The dashed lines which vanish at = 1.69 atm 
( 2 ) 

and Per = 0.11 atm correspond to the inversion frequencies 
of pure species. The critical pressure of the mixture, point 
where hoJ- vanishes, is P cr = 0.70 atm. Energy is measured 
in cm -1 . 



P (atm) 


FIG. 6 : (Color online) As in Fig. 5 for x 2 = 0.0001. For 
such small value of x 2 the critical pressure of the mixture is 
P cr ~ Per' 1 and we have hu 7+ ~ A Ei_. Inset is a zoomed in 
view of the region P ~ P cr . 
























11 


Neglecting deexcitation from these modes, i.e., sponta¬ 
neous and stimulated reemission of radiation, the absorp¬ 
tion coefficient is 

= ^ Im x(w), (38) 

with c speed of light and x( w ) as i n Eq. (32). The imagi¬ 
nary part of susceptibility is evaluated with the prescrip¬ 
tion ui —> u — irj where rj —>• 0 + . The result is 


k(lj) =1+ [S(Huj — Hu> + ) + S(Huj + fru>+)\ (39) 

+ I- [5{hjj — hoJ-) + S(Hoj + 7s*J_)], 


where the intensities I± are 

= 7 T p C(lu±) 

± 2eo ch h 2 u\. — h 2 uj‘^ z ’ 


(40) 


with C(oj) given by Eq. (33). These intensities depend 
on the pressure and for P —> 0 at the leading order they 
are 


I±{P) = x i± 


nAE i± Hj ± ^ 
Eochk bT 


(41) 


We have assumed, as usual, p = P/(k B T) and used the 
property fujJ± —> A Ei ± for P —>■ 0. Thus, when the 
density vanishes and the inter-molecular interactions can 
be neglected, the absorption coefficient of the mixture is 
the sum of the absorption coefficients of the two isolated 
species weighted by the corresponding fractions. 

Note that for non-interacting molecules, i.e., 7 ij = 0, 
i,j = 1,2, the intensities I± are as in Eq. (41), i.e., the 
absorption coefficient grows indefinitely with the pres¬ 
sure. On the other hand, for interacting molecules we 
obtain that for P —> 00 the absorption coefficient ap¬ 
proaches a constant value with asymptotic intensities 


_ x ^A^p, 2 p(1) 711 


/+ = X! 


Eochk B T 
+ x 2 


nAE 2 pl p(2) 


711^1 + 712^2 
722 


1 - 


( ae 2 


1 - 


£ 0 chk B T ~ 7 22 x 2 -I- 721 X 1 

2 (A E- 


- x, nAEWl P (1) 711 


\huj ® 
2 


E 0 cHk B T cr 711X1+712X2 \Hu 


+ 


+ x 2 


nAE 2 p2 


p{ 2 ) 


722 


/ AEi 


£ 0 chk B T cr 722X2+721X1 


f AEi x 2 


To understand the limit of the absorption coefficient 
for, e.g., x 2 -+ 0 , observe that when the fraction of species 
2 vanishes we have P cr —> PcP and the limit frequencies 
are 


Taxi .4 


hoj- 


AEiyJl- P/P£\ 
A E 2 , 

A E 2 , P < P c ( r 1} 

0, P > P c ( r 1} 


P < Pc ( r } 

P > Pc ( r } 



P (atm) 


FIG. 7: (Color online) Intensities I± of the absorption coef¬ 
ficient as a function of pressure P for the mixture of Fig. 5. 
Energy is measured in cm -1 and intensities, which have di¬ 
mensions energy divided by length, are given in cm -2 . 


if AEi > AE 2 , whereas for AE- 2 

> AEi we have 

Hoj + —► A e 2 , 




~ P / Per ' 

P < Pc^ 

10 , 


P > Pc ( r 1} 

From Eq. (40) it follows 

that 


j . ttAPiMi J 

' P, 

P < P£ ] 

EochksT | 

,°> 

P>Pc ( r } ’ 

j v ttAPiMi I 

f 0, 

P < Pc ( r 1} 

EochksT I 

[ Pc r\ 

P>Pc ( r 1} ’ 

if APi > AE 2 , whereas for AE 2 

> AUi we have 

I + -+ 0, 



T . x-AEipl 

( P 

P < Pcr^ 

EochksT 

{ P£ } 

P>Pg ) ' 


In both cases, we recover the intensity of the pure chiral 
gas 1, increasing linearly with P in the delocalized phase 
and locked at a constant value in the localized one. 

In the limit Xi —> 0 one has similar formulas which 
can be obtained from the above ones by exchanging the 
indices 1 «+ 2 . 

An example of the behavior of the intensities I± as 
a function of pressure P is illustrated in Figs. 7 and 8 
for the same NH 3 -ND 3 mixtures considered in Figs. 5 
and 6 , respectively. In this case we have AFi > A E 2 
and the fraction x 2 decreases form 0.1, Figs. 5 and 7, to 
0.0001, Figs. 6 and 8 . The dominant intensity, I + or 
is mainly due to the more abundant species i = 1 . 

We now discuss the limit in which AE\,AE 2 -+ A E, 
Mi, M 2 —1► M an d 7 i,j -+ 7 , i,j = 1,2, whereas Xi and x 2 
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P (atm) 

FIG. 8 : (Color online) As in Fig. 7 for the mixture of Fig. 6 . 


are arbitrary with X\ + X 2 = 1. In this case we have 
hu> + — » AE, 
huj- 


f AE^Jl - P/P CI , 
10 , 


P<P C , 

P>P C1 


where P cr = AE/(2 r y). It is simple to verify that 


I + — > 0 , 

nAEfj 2 I p, P< P CI 

£ 0 chk B T \ Per, P > Per ’ 

which is again the result for a single species characterized 
by the parameters A E, /r and 7 . 

Due to line broadening, which in the microwave re¬ 
gion of interest for a NH 3 -ND 3 mixture is mostly due to 
molecular scattering, the absorption coefficient discussed 
above can not be directly compared with that measured 
in an experiment. However, a test of the excitation ener¬ 
gies (36) could be done in an indirect way as in [9-11]. In 
these experiments the spectrum of the absorption coeffi¬ 
cient measured for NH 3 or ND 3 at different pressures was 
fitted by a Van Vleck-Weisskopf formula [23] in which 
the center line frequency is considered a fitting parame¬ 
ter. The results of the fit are in excellent agreement with 
the predictions of our model for a single chiral gas, see [7] 
or Appendix C. In the case of a NH 3 -ND 3 mixture the 
situation is not different. If our model is correct, with 
a 10% ND 3 fraction one should measure an absorption 
spectrum, mainly due to the hoJ+ excitation mode, see 
Fig. 7, centered at a frequency which at high pressures 
approaches the value of 0.26 cm -1 , see Fig. 5. This is 
a striking difference with respect to the case of a pure 
NH 3 gas in which the center line frequency of the ab¬ 
sorption spectrum vanishes for P > 1.69 atm. Also the 
behavior of fiw_(P) is rather different from the inversion 
frequency of pure ND 3 . In fact, the two frequencies van¬ 
ish at 0.70 atm and 0.11 atm, respectively. However, this 
effect should be more difficult to measure as the intensity 
/_ is much lower than /+, see Fig. 7. 


VI. FINAL REMARKS 


Our approach to the existence of chiral molecules is 
based on ideas of equilibrium statistical mechanics. One 
may be surprised by the presence of a quantum phase 
transition at room temperatures. We emphasize that the 
transition takes place only in the inversion degrees of 
freedom. The dynamics of these degrees of freedom is 
affected by temperature only through the values of the 
coupling constants. A derivation of our approach from 
first principles is an open question. As we have discussed, 
in the case of a pure gas of ammonia or deuterated am¬ 
monia, our model compares very well with experimental 
results. 

The problem of existence of chiral molecules has been 
approached also from a dynamical point of view [24-26] 
in terms of the quantum linear evolution with a Lindblad 
term. The question asked by these authors is how long a 
molecule in a chiral localized state remains in this state 
considering the interaction with the environment. This 
question had already been posed by Hund [1] without, 
however, any quantitative estimate. In Refs. [24, 25], a 
calculation ab initio is discussed exploiting a mechanism 
whose basic idea goes back to [27, 28]: through repeated 
scatterings in a host gas, a molecule is blocked in a lo¬ 
calized state. A detailed analysis of the decoherence of 
a two-state tunneling molecule due to collisions with a 
host gas is given in [26] in terms of a succession of quan¬ 
tum states of the molecule satisfying the conditions for a 
consistent family of histories. 

Effective nonlinearities, closer to our interpretation, 
were considered in a dynamical context by [29] and qual¬ 
itatively connected to the disappearance of the inversion 
line. However, no quantitative estimate was attempted. 
In [30] dissipative effects are included in a nonlinear time- 
dependent model obtaining qualitative conclusions simi¬ 
lar to [7]. What is missing is a theory of the parameters 
of the model. 

An important test of the present theory is to verify the 
shift of the critical pressure depending on the percentage 
of the two species in a binary mixture, e.g., NH 3 -ND 3 . 
As we have shown an addition of 10% of ND 3 is sufficient 
to halve the critical pressure of pure NH 3 and to appre¬ 
ciably change the dependence on pressure of the inversion 
frequencies of both species. 

The interpretation of the experimental data [9-11] is 
based on the Van Vleck-Weisskopf formula. The inver¬ 
sion line measured at sufficiently high pressures consists 
actually of a band of lines close to each other arising 
from different angular momentum states. At low pres¬ 
sures these lines can be resolved [31]. We may ask the 
question whether this case can be discussed in terms of a 
mixture of different pre-chiral gases corresponding to am¬ 
monia molecules in different angular momentum states 
and study the crossover to higher pressures. 
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Appendix A: Stationary states of a chiral gas 


We call these states chiral in the sense that ip\ 3 = cU-i/bu 
and i[>\ 4 = <J x ipx 3 ■ In the limit G A E, the eigenstates 
ip\ 3 and i/>x 4 approach the localized states x/Nip^ and 
y/NipR, respectively. 

The states ip\ determined above, are the stationary 
solutions ip(t) = exp(—i\t/h)ip\ of the time-dependent 
nonlinear Schrodinger equation 


The nonlinear eigenvalue problem associated with 
Eq. (2) and the normalization condition, namely, 

h(4>\)4>\ = \ip\, (VwV’a) = N, (Al) 

where 

hW = ~a x -^,a^)a z , 

has different solutions depending on the value of the ratio 
G/AE. To see it, let us write 

ipx = axy+ + bx (fi- ■ 

The coefficients ax and bx can be chosen real and are 
normalized to the number of molecules in the gas, a\ + 
b\ = N. By inserting this expression into Eq. (Al), we 
find that ax, bx are the solutions of the equation 

AEaxbx = ^ (a\ - b\) 2a x bx- (A2) 

Once a solution of this system is found, the corresponding 
eigenvalue is given by 

AM =-2~ K - h \ ) ^ (/7 a A & A- 

For any value of G, Eq. (A2) always admits solutions 
such that axbx = 0. Up to an irrelevant sign, there are 
two different solutions which satisfy this condition and 
the normalization rule, namely, the delocalized eigen¬ 
states 


Vbu = y/N<p+, Ai = -AE/2, (A3) 

V>a 2 = VN<p-, A 2 = +AE/2. (A4) 

If G > AE/2, there appear two further solutions such 
that axbx ^ 0, namely, the roots of 

2a\/N — 1 = AE/(2G). (A5) 

The corresponding eigenstates are 



A 4 = ~G, (A7) 


and the eigenvalue A is the Lagrange multiplier fixing the 
conservation of the number of molecules. Associated with 
any solution of this equation there is a conserved 
energy given by 

A 77 

The value of this functional calculated at the stationary 
solutions (A3)-(A7), i.e., 

c( I \ (2 i 2 \ C 2G 2 2 

= - jf- ( a A ^ b x) - 


provides the corresponding single-molecule energies e* = 

ety\i)/N, 4, 


AE 


AE 

e 2 =+ — 

AE 

e 3 =e 4 = - — 



These energies are plotted in Fig. 9 as a function of G. 
The molecular state effectively assumed by the gas is that 
with the minimal energy, namely, the symmetric delocal¬ 
ized state Vai for G < AE/2, or one of the two degener¬ 
ate chiral states 'ip\ 3 , ipx A for G > AE/2. 

The above results imply a bifurcation of the mean- 
field ground state at a critical interaction strength G cr = 
AE/2. According to Eq. (3), this transition can be ob¬ 
tained for a given molecular species by increasing the gas 
pressure above the critical value P cr given by Eq. (4). 


Appendix B: Excitations from the ground state 

It is interesting to compare the mean-field molecular 
states described in the previous section with the eigen¬ 
states of the A-body model (1). In the latter case, the 
Hamiltonian is represented by a 2 N x 2 N matrix which, 
if N is not too large, can be diagonalized numerically. 
For a meaningful comparison with the mean field we set 
Qij = G/N. The eigenvalues obtained for TV = 13 are 
shown in Fig. 9 as a function of G. We see that there is 
a strict correspondence between the smallest and largest 
eigenvalues of (1) and the smallest and largest conserved 
energies associated with the mean-field eigenstates of (2). 
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FIG. 10: (Color online) Energy gap E 2 — E\ of the JV-body 
Hamiltonian (1) with Qij = G/N and N = 3, 5, 7, 9, 11, 
13 (dashed lines from top to bottom) as a function of the 
ratio G/AE. The solid line represents hu/AE, where v is the 
inversion doubling frequency predicted by Eq. (B2). 


FIG. 9: (Color online) Eigenvalues E^, k = 1,...,2 N , di¬ 
vided by N, of the IV-body Hamiltonian of Eq. (1) with 
jPj = G/N and N = 13 calculated numerically for differ¬ 
ent values of G (dots). The solid lines are the single-molecule 
energies a = £(ipXi)/N associated with the mean-held sta¬ 
tionary states ipXi, i — 1,..., 4. 


Assuming that the eigenstates of (1) are ordered accord¬ 
ing to their magnitude, E\ < E 2 < ■ ■., E 2 n , we have 

/ £(V>aJ, G < AE/2, 

1 ~ 1 £ W>a 3 ,a 4 ), G > AE/2, 

E 2 n ~ £(ip a 2 ). 

By increasing the size N of the system we find that the 
differences between the IV-body and mean-field energies 
indicated above increase slower than N. Assuming that 
these differences are o(N), in Fig. 9 we compare single¬ 
particle rescaled values, namely, Ek/N, k = 1, ...,2 JV , 
and e* = £{ipx J/N, i = 1,...,4. The results shown 
for N = 13 already give an idea of what to expect for 
N —1 00 . 

The diagonalization of the Hamiltonian (1) provides 
a full description of the excitation modes. The energy 
gap between the ground state and the first excited one, 
namely, E 2 — E\, is the inversion doubling, namely, the 
excitation of one quasi-molecule from the state to 
the state </?_. For G = 0, the gap amounts exactly to 
A E and the excited level is iV-fold degenerate. Due to 
the attractive dipole-dipole interaction, the value of the 
gap evaluated for G > 0 is decreased with respect to 
the noninteracting case by an amount of the order of G. 
For G large the ground state becomes two-fold degenerate 
and the gap vanishes. In Fig. 10 we show the value of E 2 — 
E 1 evaluated as a function of G for systems of increasing 
size N. For any value of N the gap vanishes smoothly 


at G > A E, however, by increasing N such smoothness 
becomes paler and paler. What is the behavior that we 
should expect in the N —> 00 limit? As we will see in a 
moment, the mean-field indicates that E 2 — E\ sharply 
vanishes at G = AE/2 thus suggesting that for N —>■ 00 
a phase transition takes place. 

In general, the excitation energies from the ground 
state can be evaluated as the poles of the Fourier trans¬ 
form of the linear-response function [21]. In [7], we 
have evaluated the linear response of the mean-field 
Hamiltonian (2) to a time-dependent perturbation which 
corresponds to an electro-magnetic radiation coupling 
with the electric dipole of the molecules. For a gas in 
the mean-field ip \ 1 , the Fourier transform of the linear- 
response function is [32] 

2/±e 

= (M 2 - (A# 2 - 2GAE) ’ 

which has two simple poles at angular frequency zbcJ, 
id = a/A E 2 — 2GAE/H. We conclude that, in the region 
0 < G < AE/2 where the gas is at equilibrium in the 
molecular state ip\ 1 , the mean-field model predicts an 
inversion doubling frequency V = uj/(2n) given by 


v = 


A E f 2G 

"TTV a e' 


(B2) 


which vanishes at G = AE/2. This behavior of hv is 
compatible with the N —> 00 limit of the gap E 2 — E\ 
which can be envisaged from Fig. 10. 


Appendix C: Pressure dependence of NH 3 and ND 3 
inversion lines: comparison with experiments 

In [7], we have compared the mean-field theoretical 
prediction for the inversion frequency with the spectro- 
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FIG. 11: (Color online) Plot of the universal function of 
Eq. (Cl). Points are the experimental results for NH 3 (cir¬ 
cles) or ND 3 (squares) available from Ref. [9, 10] and [11], 
respectively. The experimental data have been scaled using 
the values of A E and P CI from Table I. 


scopic data available for ammonia [9, 10] and deuterated 
ammonia [11]. In these experiments, the central fre¬ 
quency of the inversion peak was measured at different 
pressures. The frequency is found to decrease by increas¬ 
ing P and vanishes for pressures greater than a critical 


value. This behavior is very well accounted for by the 
mean-field prediction (B2) for v. In fact, an important 
feature of this equation is the statement that the ratio 
between the frequency V of a measured spectroscopic in¬ 
version line and the energy splitting A E of an isolated 
molecule, 


— = h- — 

A E V Per’ 


(Cl) 


is a universal function of P/P cr independent of the kind of 
chiral species considered. In Fig. 11, we show the values 
of hv/ AE measured at different pressures P in the case 
of ammonia and deuterated ammonia. The superposition 
of the two sets of data plotted as a function of P/P CI and 
the agreement with formula (Cl) is impressive. 

Some comments are in order. The mean-field predic¬ 
tion (Cl) does not contain free parameters. The value of 
P cr is evaluated from the electric dipole p of the molecule, 
its collision diameter d and inversion energy splitting A E, 
and the temperature T of the gas, see Eq. (4) and Table I. 
We remark, however, that the experimental points shown 
in Fig. 11 are the result of an indirect measurement. 
More precisely, in experiments [9-11] the line shape of 
the absorption coefficient measured for NH 3 and ND 3 at 
different pressures P is fitted by a Van Vleck-Weisskopf 
formula [23] in which V is a considered a fitting parame¬ 
ter. 
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